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Abstract: The theoretical description of nuclear resonances at zero and finite temperatures 
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CN ■ is presented. The following issues are addressed: 

1) Giant dipole resonances (GDR) in highly excited nuclei, including both low and high 
regions of temperature. The results of calculations are obtained within the phonon-damping 
model, thermal shape-fluctuation model including thermal pairing, and compared with ex- 
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^ \ perimental data. 

2) The electromagnetic cross sections of the double GDRs (DGDR) in 136 Xe and 208 Pb. 
The results obtained in theoretical calculations are compared with the experimental data 
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for the DGDR cross sections in exclusive measurements at near-relativistic energies. 

3) GDR and pygmy dipole resonances (PDR) in neutron-rich nuclei, where the effect of 
coupling of the GDR to complicated configurations on the PDR is analyzed. 

1. Introduction 

Nuclear resonances under extreme conditions mentioned in the present lecture include 
i) giant dipole resonances (GDR) in highly excited nuclei formed in heavy-ion reactions and 
inelastic scattering of light particles or nuclei on heavy targets, ii) double GDRs (DGDR) 
formed in Coulomb excitations at near relativistic energies, and iii) pygmy dipole resonances 
(PDR) in neutron-rich nuclei. 

The recent studies of these resonances have been facing the following challenges: 
1) The GDR has been observed in highly-excited (hot) nuclei. These nuclei are formed 
as compound nuclei at high excitation energies E* in heavy-ion fusion reactions or in the 
inelastic scattering of light particles (nuclei) on a heavy nucleus. The 7-decay spectra of 
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these compound nuclei show the existence of the GDR, whose peak energy depends weakly 
on the excitation energy E*. The dependence of the GDR on the temperature T has been 
experimentally extracted when the angular momentum of the compound nucleus is low, as 
in the case of the light-particle scattering experiments, or when it can be separated out from 
the excitation energy E*. It has been experimentally found that the GDR's full width at 
the half maximum (FWHM) remains almost constant at T < 1 MeV, but sharply increases 
with T up to T ~ 2 - 3 MeV, and saturates at higher T 

2) The DGDR has been observed in the Coulomb excitations of 136 Xe and 208 Pb projectiles 
on 208 Pb target at the bombarding energies of 0.7 GeV/n and 0.64 GeV/n, respectively. 
The extracted electromagnetic (EM) cross sections of the DGDR in these experiments have 
been compared with theoretical predictions based on the non-interacting phonon (harmonic) 
picture. The latter calculates the parameters of the DGDR by folding to independent GDRs. 
As a result, the DGDR energy ^dgdr is found to be equal to 2E GDR (-^gdr is the GDR 
energy), and the DGDR FWHM r DGD R is equal to 2r GD R (r GD R is the GDR FWHM) if 
Lorentzian photoabsorption cross sections are used in folding the GDRs, or to V^Fqdr, if 
Gaussians are folded. The comparison shows that the experimentally extracted energy and 
width of DGDR differ slightly from these values. The most striking discrepancy is that the 
experimental EM cross sections are much larger that than the values predicted by the folding 
model. The enhancement is found to be around 178% for 136 Xe and 133% for 208 Pb j^J. 

3) It is well believed that, in neutron-rich nuclei, the oscillation of the neutron excess 
against the stable core may cause a low-frequency El resonance, which is called the pygmy 
dipole resonance (PDR). However, in the case when the neutron excess cannot be well sep- 
arated from the stable core, the coupling between them may lead to a strong fragmentation 
of the PDR. As a result the pygmy dipole mode may become non-collective. Theoretical cal- 
culations within the relativistic RPA and the particle- hole (p/i)©phonon-coupling model so 
far have given contradicting predictions. The former shows a clear PDR bump with a peak 
at around 8 MeV, while no evidence for collectivity is seen below 10 MeV in the latter jjjj]. 

It will be shown in the present lecture how the above-mentioned issues are interpreted 
within the phonon-damping model (PDM) The PDM was originally proposed in 1998 
with the primary aim to describe the damping of the hot GDR. This model was extended 
later to address the issue of anharmonicity in the DGDR. By including pairing, the model 
is also able to predict the GDR's width at low T as well as the EM cross sections of El 
excitations in neutron-rich nuclei. 
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2. Outline of the PDM 



A. Damping of hot GDR 

The quasiparticle representation of the PDM Hamiltonian \^ is obtained by adding the 
superfluid pairing interaction and expressing the particle (p) and hole (h) creation and 
destruction operators, a\ and a s (s — p,h), in terms of the quasiparticle operators, ot\ and 
a s , using the Bogolyubov's canonical transformation. As a result, the PDM Hamiltonian 
for the description of EA excitations can be written in spherical basis as 
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where A = a/2A + 1. The first term at the right-hand side (rhs) of Hamiltonian ([!} cor- 
responds to the independent-quasiparticle field. The second term stands for the phonon 
field described by phonon operators, b\^ and b\^i, with multipolarity A, which generate the 
harmonic collective vibrations such as GDR. Phonons are ideal bosons within the PDM, i.e. 
they have no fermion structure. The last term is the coupling between quasiparticle and 
phonon fields, which is responsible for the microscopic damping of collective excitations. 
In Eq. ((H) the following standard notations are used 
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combinations of Bogolyubov's u and v coefficients. The quasiparticle energy Ej is calculated 
from the single-particle energy 6j as 
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where the pairing gap A and the Fermi energy e F are defined as solutions of the BCS 
equations. At T ^ the thermal pairing gap A(T) (or A(T)) is defined from the finite- 
temperature BCS (or modified BCS) equations (See section C below). 



The equation for the propagation of the GDR phonon, which is damped due to coupling 
to the quasiparticle field, is derived making use of the double-time Green's function method 
(introduced by Bogolyubov and Tyablikov, and developed further by Zubarev). Follow- 
ing the standard procedure of deriving the equation for the double-time retarded Green's 
function with respect to the Hamiltonian ((TJ, one obtains a closed set of equations for the 
Green's functions for phonon and quasiparticle propagators. Making the Fourier transform 
into the energy plane E, and expressing all the Green functions in the set in terms of the 
one-phonon propagation Green function, we obtain the equation for the latter, G\i(E), in 
the form 

where the explicit form of the polarization operator P\i(E) is 
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The polarization operator (jSJ) appears due to ph - phonon coupling in the last term of the 
rhs of Hamiltonian The phonon damping 7Ai(w) (uj real) is obtained as the imaginary 
part of the analytic continuation of the polarization operator P\i(E) into the complex energy 
plane E = uo ± is. Its final form is 
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{v { r)f{ nj - n f )[5(E - E j + E r ) - 8{E + e 3 - e,-,)]}- 
The quasiparticle occupation number rij is calculated as 

n J-oo [E - Ej - Mj(E)] 2 + Yj{E) 

where Mj(E) is the mass operator, and ^fj{E) is the quasiparticle damping, which is de- 
termined as the imaginary part of the complex continuation of Mj(E) into the complex 
energy plane [5]. These quantities appear due to coupling between quasiparticles and the 
GDR. From Eq. (|7J) it is seen that the functional form for the occupation number n - is not 
given by the Fermi-Dirac distribution rip(Ej) for non-interacting quasiparticles. It can be 
approximately to be so if the quasiparticle damping jj(E) is sufficiently small. Equation ((7|) 
also implies a zero value for rij in the ground state, i.e. rij{T = 0) = 0. In general, it is not 
the case because of ground-state correlations beyond the quasiparticle RPA (QRPA). They 



lead to rij(T = 0) 7^ 0, which should be found by solving self-consistently a set of nonlinear 
equations within the renormalized QRPA. However, for collective high-lying excitations such 
as GDR, the value rij(T = 0) is negligible. 

The energy uj of giant resonance (damped collective phonon) is found as the solution of 
the equation: u — u\i — P\i(w) = . The width T\ of giant resonance is calculated as twice 
of the damping 7 A (u) at uj = Co, where A = 1 corresponds to the GDR width Tgdr- The 
latter is conveniently decomposed into the quantal (Fq) and thermal (r T ) widths as 

r G DR = r Q + r T , (8a) 
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where (ss') = (pp 1 ) and (hh') with p and h denoting the orbital angular momenta j p and 
jh for particles and holes, respectively. The quantal and thermal widths come from the 
couplings of quasiparticle pairs [ad <8> «[]im and [aj <8> a-^LM to the GDR, respectively. At 
zero pairing they correspond to the couplings of ph pairs, [al ® aj^LM, and pp (hh) pairs, 
[a\ ® a~,]LM, to the GDR, respectively (The tilde ~ denotes the time-reversal operation). 

The line shape of the GDR is described by the strength function SgdrO^), which is 
derived from the spectral intensity in the standard way using the analytic continuation of 
the Green function (J3J) and by expanding the polarization operator (JHJ) around u = -Eqdr- 
The final form of Sgdr(^) is j^, |5j] 
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The photoabsorption cross section cr(£' 7 ) is calculated from the strength function S'gdrC^) 
as 

a(E 7 ) = Ci> Sgdr(^ 7 )^ 7 , (10) 

where E 1 = u is used to denote the energy of 7-emission. The normalization factor C\ is 
defined so that the total integrated photoabsorption cross section a — J cr(E 1 )dE 1 satisfies 
the GDR sum rule SRgdr, hence 
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c\ = SRgdr / / SGr>n(E 7 )E 7 dE 7 . (11) 



In heavy nuclei with A > 40, the GDR exhausts the Thomas- Reich-Kuhn sum rule (TRK) 
SRqdr = TRK = 60 NZ/A (MeV-mb) at the upper integration limit £ max ~ 30 MeV, and 
exceeds TRK (SRqdr > TRK) at E iaax > 30 MeV due to the contribution of exchange 
forces. In some light nuclei, such as 16 0, the observed photoabsorption cross section 
exhausts only around 60% of TRK up to E max ~ 30 MeV. 

B. Thermal pairing 

The standard finite-temperature BCS (FT-BCS) theory ignores fluctuations of the quasi- 
particle number. As a result, the BCS breaks down at a critical temperature T c ~ 
0.567A(T = 0), which corresponds to the sharp transition from the superfluid phase to 
the normal-fluid one. It has been known that, in finite systems such as nuclei, thermal 
fluctuations smooth out this phase transition 61. 

ri 

The modified BCS (MBCS) theory 7] proposes a microscopic way to include quasiparticle- 
number fluctuations via the secondary Bogolyubov's transformation 

«jm = a/1 - njOi] m + yfn~ajm , a jfh = ^/n~a jfh - y/n~j(x] m . (12) 

Using Eqs. (fT2"|) in combination with the original Bogolyubov's transformation, one ob- 
tains the transformation from the particle operators directly to the modified quasiparticle 
operators in the following form 

a ]m = Uj&jm + VjOtjm , (Ijfc = UjOljf^ — VjOij m , (13) 

where the coefficients Uj and Vj are related to the conventional Bogolyubov's coefficients Uj 
and Vj as 

Uj = Ujy/l — rij + Vjy/rTj , Vj = Vj^Jl — rij — Uj^/rTj . (14) 

The transformation of the pairing Hamiltonian (|1J) into the modified quasiparticles aj m and 
OLj m has the form identical to that obtained within the conventional quasiparticle represen- 
tation with (Uj, Vj) replacing (uj, Vj) and (aj m , atj m ) replacing (a] m , <x, m ), respectively. 
The MBCS equations, therefore, has exactly the same form as that of the standard BCS 
equations, where the coefficients Uj and Vj are replaced with Uj and Vj, i.e. 

A = G ttjUjVj = Gj2 n A( 1 - 2n j) u J v j ~ \l n ji l - n i) K 2 - v ])\ > ( 15 ) 
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The last terms at the rhs of these MBCS equations contain the quasiparticle-number fluctu- 
ations y/rij(l — rij) on j— th orbitals, which are not included in the standard FT-BCS theory. 



C. EM cross sections of GDR and DGDR 

The EM cross section oem is calculated from the corresponding photoabsorption cross 
section cr(E 7 ) and the photon spectral function N(E 7 ) as 

/poo 
N{E 1 )a{E 1 )dE 11 iV(£ 7 ) = 2vr / e" m(b) N(E^,b)bdb. (17) 
•'f'min 

The expression for the spectrum N(E 7 , b) of virtual photons from a stationary target as 
seen by aprojectile moving with a velocity f3 = v/c at the impact parameter b is also 
given in |8(. The average number of photons absorbed by the projectile is calculated as 
m{b) = J£ ;n iV(£ 7 ,&)<x(£ 7 )d£ 7 . 

The DGDR strength function is calculated within the PDM as 

qPdm /en 2 7dgdr(£) /1cA 
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where the DGDR energy -Edgdr and damping 7dgdr(-E) are calculated microscopically 
within the PDM (See details in Ref. [9|). The DGDR cross section ctdgdr^) is calculated 

as 

odgdr(-E') = c®S™ R (E)E . (19) 

The DGDR strength factor c?' in (|19|) as follows. Using ()17jl and the harmonic limit 
^DGDR(har) (-^) °f the DGDR strength function (JTHJ), which is obtained by folding two GDR 
strength functions (pairing not included), we write the formal expression of the harmonic 
limit a^(har) of the EM cross section (JT7J) for DGDR as 

/A (^) P 
^(har)rfE = J NUE)S£™ R{hai) (E)EdE, (20) 

where N har (E) is calculated using the harmonic limit fDGDR(har) (E) of (|T9"j) in (fTT|) and m{b). 
We require this cross section (ffil)) to be equal to the one calculated by folding two GDR 
cross sections, namely 
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where the representation 5(x) = [(x — is)^ 1 — (x + ie)~ 1 ]/(2'iri) and the expression for 
N(Ei, E 2 ) given in [8j are used. Equalizing the right-hand sides of (f2T)|) and (|2*T|) . we define 
c^ 2 \ Knowing c^ 2 \ we can calculate the EM cross section of the DGDR from (J2*U|) using 
^dgdr(^) ffTH|) instead of its harmonic limit. 

3. Numerical results 

A. Assumptions and parameters of PDM 

The PDM is based on the following assumptions: 

al) The matrix elements for the coupling of GDR to non- collective ph configurations, 
which causes the quantal width Tq (|8b[). are all equal to Fx. Those for the coupling of GDR 
to pp (hh), which causes the thermal width i\ (jHcjl . are all equal to F 2 . 

a2) It is well established that the microscopic mechanism of the quantal (spreading) width 
Tq (|8b[) comes from quantal coupling of ph configurations to more complicated ones, such 
as 2p2h ones. The calculations performed in Refs. within two independent microscopic 
models, where such couplings to 2p2h configurations were explicitly included, have shown 
that Tq depends weakly on T. Therefore, in order to avoid complicate numerical calculations, 
which are not essential for the increase of Tgdr at T ^ 0, such microscopic mechanism is 
not included within PDM, assuming that Tq at T = is known. The model parameters are 
then chosen so that the calculated Tq and -Egdr reproduce the corresponding experimental 
values at T = 0. 

Within assumptions (al) and (a2) the model has only three T-independent parameters, 
which are the unperturbed phonon energy u q , Fx, and F 2 . The parameters u q and Fx are 
chosen so that after the p/i-GDR coupling is switched on, the calculated GDR energy -Egdr 
and width Tgdr reproduce the corresponding experimental values for GDR on the ground- 
state. At T ^ 0, the coupling to pp and hh configurations is activated. The F 2 parameter 
is then fixed at T = so that the GDR energy -Eqdr does not change appreciably with T. 

B. Temperature dependence of GDR width 

Shown in Fig. Q(a) is the T dependence of the neutron pairing gap A v for 120 Sn, which is 
obtained from the MBCS equation (fTKj) using the single-particle energies determined within 
the Woods-Saxon potential at T = 0. The pairing parameter G v is chosen to be equal to 
0.13 MeV, which yields A(T = 0) = A(0) ~ 1.4 MeV. Contrary to the BCS gap (dotted 
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FIG. 1: (a): Neutron pairing gap for 120 Sn as a function of T. Solid and dotted lines show the 
MBCS and BCS gaps, respectively, (b): GDR width Tgdr as a function of T for 120 Sn. The 
thin and thick solid lines show the PDM results obtained neglecting pairing and including the 
renormalized gap A = [1 + 1/5N 2 ]A, respectively. The gap A includes the correction 5N 2 = 
A(0) 2 + V^)/[( e j ~~ ^f) 2 + A(0) 2 ] due to an approximate number projection. The prediction 
by the TFM is shown as the dotted line 

line), which collapses at T c ~ 0.79 MeV, the gap A (solid line) does not vanish, but decreases 
monotonously with increasing T at T > 1 MeV resulting in a long tail up to T ~ 5 MeV. 
This behavior is caused by the thermal fluctuation of quasiparticle number in the MBCS 
equations (|T5|l and (JTBJ). 

The GDR widths as a function of T for 120 Sn obtained within the PDM are compared in 
Fig. (b) with the experimental data and the prediction by the thermal fluctuation model 
(TFM) [11]. The TFM interprets the broadening of the GDR's width via an adiabatic 
coupling of GDR to quadrupole deformations induced by thermal fluctuations. Even when 
thermal pairing is neglected the PDM prediction, (the thin solid line) is already better than 
that given by the TFM, including the region of high T where the width's saturation is 
reported. The increase of the total width with T is driven by the increase of the thermal 
width Tt (JHcj) . which is caused by coupling to pp and hh configurations, since the quantal 
width Tq (J8bj) is found to decrease slightly with increasing T. The inclusion of thermal 
pairing, which yields a sharper Fermi surface, compensates the smoothing of the Fermi 
surface with increasing T. This leads to a much weaker T-dependence of the GDR's width at 
low T. As a result, the values of the width predicted by the PDM in this region significantly 
drop (the thick solid line), recovering the data point at T = 1 MeV. 

The results discussed above have also been confirmed by our recent calculations within a 
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FIG. 2: (a): Pairing gaps for 120 Sn averaged over thermal shape fluctuations versus T. Lines with 
triangles and crosses are the usual BCS proton and neutron pairing gaps, respectively, while those 
with diamonds and squares denote the corresponding pairing gaps, which also include thermal 
fluctuations of pairing fields, (b): GDR widths for 120 Sn versus T. Open squares, triangles, 
and diamonds denote the widths obtained without pairing, including BCS pairing, and thermally 
fluctuating pairing field from (a), respectively. 

macroscopic approach, which takes pairing fluctuations into account along with the thermal 
shape fluctuations [12]. Here the free energies are calculated using the Nilsson-Strutinsky 
method at T ^ 0, including thermal pairing correlations. The GDR is coupled to the nuclear 
shapes through a simple anisotropic harmonic oscillator model with a separable dipole-dipole 
interaction. The observables are averaged over the shape parameters and pairing gap. Our 
study reveals that the observed quenching of GDR width at low T in 120 Sn and 148 Au can 
be understood in terms of simple shape effects caused by pairing correlations. Fluctuations 
in pairing field lead to a slowly vanishing pairing gap [Fig. [21 (a)], which influences the 
structural properties even at moderate T (~1 MeV). We found that the low-T structure and 
hence the GDR width are quite sensitive to the change of the pairing field [Fig. [21 (b)]. 
C. EM cross section of DGDR 

The calculations employ the single-particle energies for 136 Xe and 208 Pb obtained within 
the Hartree-Fock (HF) method using the SGII interaction. The unperturbed phonon energy 
uj q and the degenerate ph matrix element F p h are chosen so that the GDR energy £"gdr and 
FWHM Tgdr; obtained within the PDM, reproduce the experimentally extracted values. 

The peak energies FWHM r i; and EM cross section o~ l c for GDR (i=GDR) and 
DGDR (i=DGDR) predicted by the PDM are shown in Table [1] in comparison with the 
experimental data for 136 Xe and 208 Pb |2j. It is seen that all the calculated values are in 
reasonable agreement with the corresponding experimental values. K. Boretzky of the LAN 
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TABLE I: The energies E u FWHM r i} and EM cross section a l c for GDR (i=GDR) (a) and 
DGDR (i=DGDR) (b), calculated within PDM (Theory) in comparison with the experimental 
data (Exper.) for 136 Xe and 208 Pb 



a 


E GDR (MeV) 


r GD R (MeV) 


gdr / i \ 
<7g UK (mb) 


AT 1 

Nucleus 


Theory Exper. 


Theory Exper. 


Theory Exper. 


136 Xe 


15.6 15.2 


4.96 4.8 


1676.28 1420(42)±100 


208p b 


13.5 13.4 


4.04 4.0 


3039.67 3280±50 


b 


£dgdr (MeV) 


Tdgdr (MeV) 


^DGDR (mb) 


Nucleus 


Theory Exper. 


Theory Experiment 


Theory Exper. 


136 Xe 


29.2 28.3 ±0.7 


7.0 6.3 ±1.6 


159.33 164(85)±35 


208p^ 


26.6 26.6 ±0.8 


6.3 6.3 ±1.3 


420.92 380 ±40 
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FIG. 3: Differential EM cross section for 136 Xe (a) and 208 Pb (b). Results obtained using EM 
cross section of DGDR within PDM (solid lines) and the best fit with x 2 to the data points (dotted 
lines) are based on a normalization of GDR (the bump at ~ 10 MeV) which exhausts 90% of TRK. 



collaboration has folded the PDM strength functions for GDR and DGDR with the detector 
response and plotted the obtained results against the experimental fits in Fig. El which 
shows a remarkable agreement between the PDM predictions and the experimental data. 
D. El-excitations in neutron-rich nuclei 

The calculations of photoabsorption and EM cross sections have been carried out for 
oxygen isotopes with A = 16, 18, 20, 22, and 24, and for calcium isotopes with A = 40, 
42, 44, 46, 48, 50, 52, and 60. The calculations employ the spherical-basis single-particle 
energies Ej obtained within the HF method using the SGII interaction. The two PDM 
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E y (MeV) Ey (MeV) 

FIG. 4: Photoabsorption cross sections obtained within the PDM for some oxygen isotopes. 

parameters u>\ (A = 1) and F\ = /jj; for all ph indices (j = p,j' = h) are chosen in such 
a way that the values of GDR energy -Eqdr and width Tgdr for 16 and 40 ' 48 Ca reproduce 
their corresponding experimental values -Egdr an d r^ R . These chosen values of PDM 
parameters are then fixed in calculations for the neighbor neutron-rich isotopes (N > Z). 
The neutron pairing gap A n is adjusted around the general trend Ylj\[A of the observed 
pairing gaps in stable nuclei to keep the GDR energy stable against varying the neutron 
number N. 

Shown in Fig. 0] are the photoabsorption cross sections obtained within the PDM for 
i6,i8,2o,24Q j£ j s seen that the GDR becomes broader with increasing the neutron number 
N. Its width is particularly large for 18 ' 20 O, for which the values of neutron pairing gap A n 
are largest. This increase of GDR spreading enhances both of its low- and high-energy tails. 
In the region E 1 < 15 MeV, some weak structure of PDR is visible for 18 ' 20 O, and also 22 
(not shown). In the rest of isotopes under study, except for an extension of the GDR tail 
toward lower-energy, there is no visible structure of the PDR. 

Since the photon spectral function N(E y ) in the EM differential cross section da^yi/dE^ 
(Eq. p7j1 ) is an exponentially decreasing function with increasing E 1 , it enhances the 
low-energy part of the GDR. Therefore, the EM differential cross section da-^M/dE^ of GDR 
can be used as a magnifying glass for the structure of PDR. These cross sections obtained 
within the PDM are shown in Fig. El for 20 > 22 O. The calculations were carried out for 208 Pb 
target at various projectile energies as shown in these figures. The PDR shows up in the 
EM cross sections as a well isolated peak located at around 6 MeV. Varying the projectile 
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FIG. 5: Electromagnetic cross sections of GDR within PDM for 20 - 22 O on 208 Pb target. Different 
lines display results obtained at different projectile energies as indicated in the panels. 

energy affects mostly the GDR region (above 10 MeV) of the EM cross section, but not the 
PDR. 

4. Conclusions 

A review is given on the current status of the study of nuclear resonances. It is demon- 
strated that the PDM is a simple but yet microscopic model, which is able to describe a 
variety of resonances under extreme conditions, namely: 1) the GDR at T ^ 0, including 
the constant width at T below 1 MeV, the width increase at low T, the width saturation at 
high T; 2) the EM cross sections for DGDR obtained at near-relativistic energy in 136 Xe and 
208 Pb; 3) the possibility of extracting PDR in the EM cross sections of neutron-rich isotopes 
at low-energy beams. 

Regarding the GDR at T ^ 0, it is shown that the mechanism of the width's increase at 
1 < T < 3 MeV and its saturation at T > 3 MeV comes from the coupling of the GDR to 
non-collective pp and hh configurations at T ^ 0. Meanwhile this effect is nearly cancelled 
by the monotonous decrease of the thermal pairing with increasing T at T < 1 MeV. As the 
result, the GDR width in this low-T region remains nearly temperature independent. This 
effect is confirmed by calculations in both of the PDM and a macroscopic approach, which 
takes into account thermal fluctuations of nuclear shapes and pairing field. 

Concerning the DGDR, the PDM has succeeded to include the effect of anharmonicity 
between two coupled GDRs. As a result, for the first time, the experimental values of EM 
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cross sections are nicely reproduced by theoretical calculations within the PDM for both 
136 Xe and 208 Pb. 

Finally, the results of calculations within the PDM have demonstrated how the pygmy 
dipole mode can be depleted due to its coupling to the GDR. This can be the case when 
the neutron excess cannot be well separated from the stable core. This effect leads to the 
disappearance of collectivity of the GDR. As a result the photoabsorption cross sections 
for neutron-rich isotopes have a tail extended toward low-energy region instead of a well- 
pronounced PDR peak. Nonetheless, since the photon spectral function N(E y ) in the EM 
differential cross is an exponentially decreasing function with increasing E y , it enhances the 
low-energy part of the El-strength distribution, which is insensitive to the variation of beam 
energy. This feature suggests that a clean PDR peak (without admixture with the GDR) 
can be seen in the EM cross section (of neutron- rich oxygen and calcium isotopes, e.g.) using 
low-energy, but high-intensity beams at ~ 50 Mev/n. 
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